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Abstract 

The sensitivity of semi-hard (p± < 10 GeV) hadron production to 
parton energy loss in high energy nuclear collisions is studied via the HI- 
JING1.35 model. We test the model on recent WA98 data on 160 AGeV 
Pb + Pb — > 7T° up to 4 GeV/c and while these data are reproduced, the 
results depends sensitively on the model of the Cronin effect. At (RHIC) 
collider energies (-y/i > 200 AGeV), on the other hand, semi-hard hadron 
production becomes insensitive to the above model and is thus expected to 
be a cleaner probe of jet quenching signature associated with non-Abelian 
radiative energy loss. 

Preliminary WA98 Pb + Pb -> tt° data in the p± w 2 - 4 GeV range 
were analyzed recently by Wang Q via a parton model. The Leading Log 
Approximation (LLA) parton model |J was shown to reproduce the observed 
7T° invariant inclusive cross sections simultaneously in p + p, S + S Q , and 
central Pb + Pb |] in the CERN/SPS energy range, y/a < 20 AGeV. What 
is remarkable about that analysis is the implied absence of quenching of high 
transverse momenta hadrons that should be observed if partons lose energy in 
dense matter || ^| . The S + S and Pb + Pb data clearly reveal the expected 
Cronin enhancement of moderate p± pions, as seen first in p+A. The pp data 
also confirm the need to supplement the LLA with non-perturbative intrinsic 
fej_ (|, |^1 • However, the Pb + Pb data show no sign of jet quenching expected at 
collider energies ||. 

In central heavy ion reactions at the SPS, conservative estimates of the 
initial energy density suggest that e(l fm/c) k 1 ~ 5 GeV/fm 3 , are reached. 
The gluon radiative energy loss of partons in such dense matter is expected 
to exceed dE/dx ~ 1 GeV/fm Recent theoretical analysis JTo| predicts in 
fact a much larger non-linear energy loss AE ~ (Ax) 2 GeV/fm 2 if a parton 
traverses a quark-gluon plasma of thickness Ax. On the other hand, the WA98 
data seem to rule out dE/dx > 0.1 GeV/fm, as Wang showed in Ref. (|]. In 
this work, we consider this problem using the nuclear collision event generator, 
HIJING 



Recall that in the LLA, pQCD predicts that the single inclusive hadron cross 
section is given by Refs. [|[ || 
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This formula convolutes the elementary pQCD parton-parton cross sections, 
da(ab — > cd), with non-perturbative lowest order fits of the parton structure 
functions, f a /Ai an d the parton fragmentation functions, D h / C , to ep and e + e~ 
data. Here K a , Kb are the intrinsic transverse momenta of the colliding partons. 
The model includes a K « 2 factor to simulate next-to leading order correc- 
tions jl3) at a hard scale Q ~ pr/zc- The scale dependence of the structure 
and fragmentation functions account for the multiple soft collinear radiative ef- 
fects. However, below energies y/s < 100 GeV, it is well known m that LLA 
significantly underpredicts the p± < 10 GeV cross section, and additional non- 
perturbative effects must be introduced to bring LLA into agreement with data. 
Unfortunately, as emphasized in Ref. ||, the results are then quite model de- 
pendent below collider (RHIC) energies. As we show below, the good news is 
that this model dependence is reduced significantly at collider (RHIC) energies. 

In spite of the inherent ambiguity of the parton model analysis at SPS 
energies, a successful phenomenological approaches to this problem has been 
developed via the introduction j^] of intrinsic transverse momenta of the col- 
liding partons as in (Q). Originally, a Gaussian form for that distribution with 
~ 1 GeV 2 was proposed in Ref. Q. However, in order to reproduce pp 
data more accurately and to take into account the Cronin effect, the Gaus- 
sian ansatz was generalized in Ref. [0 to include Q 2 and A dependence. With 
<?(k) — * g a (K,Q 2 ,A), excellent fits [g] to the WA98 data could be obtained 
assuming a factorized Gaussian distribution with 

{n 2 {b))A = (l + 0.2 Q 2 a s (Q 2 ) + °- 23fT f^ ln2Q ) (GeV/c) 2 (2) 

where Q 2 is measured here in (GeV/c) 2 . In (||) a pp tA(b) is the average number of 
inelastic scatterings a nucleon suffers traversing nucleus A at impact parameter 
b. The nuclear thickness function, tA(b), is normalized as usual to J d 2 btA(b) = 
A. Eq.(||) is the main source of the model dependence in Ref. 

The HIJING1.35 Monte Carlo model |I| @ incorporates pQCD jet pro- 
duction together with initial and final state radiation according to the PYTHIA 
algorithm [H. In addition, it incorporates a variant of the soft string phe- 
nomenology similar to the Lund/FRITIOF jl5| and DPM Jl6| models to sim- 
ulate beam jet fragmentation and jet hadronization physics. Low transverse 
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momenta inelastic processes are of course highly model dependent, and the pa- 
rameters must be fit to pp and AB data O, [T^]. It is of interest to apply 
HIJING to the present study because it incorporates in addition to the above 
soft and hard dynamics, a model of soft (Cronin) multiple initial state collision 
effects as well as a simple jet quenching scheme. With these features, we are able 
to study how competing aspects of the reaction mechanism influence hadronic 
observables and explore the magnitude of theoretical uncertainties. 

In the HIJING model, excited baryon string are assumed to pick up random 
transfer momentum kicks in each inelastic scattering according to the following 
distribution 

ff ( f ?)(x{( K 2 +^)( K 2 +p 2 )(1 + e -( K ^p 2 )/p3|- 1 ; (3) 

where pi = 0.1, P2 = 1.4, pz = 0.4 GeV/c were chosen to fit low energy (p± < 1 
GeV/c) multiparticle production data in Refs. Q @. A flag, IHPR2(5)(=1 or 
0), makes it possible to compute spectra with and without this effect as shown 
in part (a) of Figs. 1-3. The present study is the first test of this model up to 4 
GeV/c in the SPS energy range. 

Jet quenching is modeled via gluon splitting according to the number of 
mean free paths, A = 1 fm, traversed by a gluon through the cylindrical nuclear 
reaction volume. In each partonic inelastic interaction a gluon of energy AE = 
Ax dE/dx is assumed to be split off the parent jet and incorporated as a kink 
in another baryonic string The (constant) energy loss per unit length is an 
input parameter (HIPR1(14) in HIJING JT1J) and can switched on and off via 
IHPR2(4) (=1 or 0) to test the sensitivity of spectra to jet quenching as shown 
in Figs. 1-3. 

Figure. |l| compares 8 the predictions of HIJING1.35 jn], [l^] without jet 
quenching (dE/dx = 0) for the invariant 7r° cross section in central nuclear 
collisions at SPS and RHIC energies. The cross section for central A + A col- 
lisions are computed integrating over the impact parameters up to & max chosen 
to reproduce experimental trigger cross sections. For WA98 Pb + Pb and RHIC 
Au + Au we took b max = 4.5 fm, while for the WA80 S + S we took & max = 3.4 
fm. The multiple collision eikonal geometry in HIJING is based on standard 
Wood-Saxon nuclear densities. 

The parton model fit to the WA98 data are labeled by 'Wang' from Ref. g . 
(The normalization of both the WA98 data and Wang's latest calculations (not 
shown) have increased ~ 20-40% relative to ||.) We note that the HIJING1.35 
calculation for this interaction given by the solid jagged line also reproduces 
remarkably well the ir° invariant cross sections without jet quenching. However, 
for the lighter S + S reaction, HIJING, underestimates the p± > 1 GeV/c tail 
significantly. This error is traced to the failure of the model to reproduce the 
pp high p± data at these energies, in contrast to its successful account of higher 
energy data|ll[]. This can be seen by comparing the filled squares to the dot- 
dashed curves as explained below. Therefore, we find that the agreement with 
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the WA98 data is accidental and the observed A scaling of the high p± region 
at SPS energies is not reproduced. 

Fig. la shows that the soft transverse momentum kick model is the source for 
the agreement of HIJING with the Pb + Pb data. The dot-dashed curves show 
what happens if the soft p±_ kicks modeled with eq. (^) is turned off. The very 
strong decrease of the pion yield in the Pb + Pb case and the somewhat smaller 
but still large decrease in the S+S case shows clearly the important role multiple 
transverse kicks at these energies. In fact both S + S and Pb + Pb dot-dashed 
curves are found to coincide with the calculated pp — > 7r differential cross section 
scaled by the wounded nucleon factor Wa&Aa/&nNi where Wa = 21, 172 is the 
average number of wounded projectile nucleons and oaa = 32,363,636 mb for 
A = 1,32,207. 

On the other hand, the data follow closely the shape of the measured pp — ► 
7r + data taken from j| and scaled to AA by multiplying by the Glauber bi- 
nary collision number factor, Taa&aa- From HIJING the average number of 
binary collisions in these systems is 45, 751 resp. The fact that the WA98 data 
scale with the above Glauber factor within a factor of three, suggests that the 
additional p± broadening due to initial state collisions is relatively small. The 
very large A dependence of the HIJING p± tail at SPS energies is due to the 
A 1 / 3 times convolution of the distribution (||). This problem is avoided in the 
parton model calculation^) using (0) through the separation of larger intrinsic 
momentum effects and smaller A 1 / dependent contributions. 

We conclude that the missing intrinsic transverse momentum component of 
HIJING precludes an accurate simultaneous reproduction of of pp, SS, PbPb 
data at SPS energies. However, unlike the parton model where no global con- 
servation laws have to be enforced, it is not clear how to incorporate intrinsic 
momenta in a global event generator like HIJING without destroying the sat- 
isfactory reproduction of low p± < 1 GeV/c data. We do not attempt to solve 
this problem here. 

Our main point is that this problem goes away fortunately at higher (RHIC) 
collider energies (y/s — 200 AGeV). As seen in Fig. la) the effect of multiple 
soft interactions is very much reduced at that energy. This is due to the well 
known effect that as the beam energy increases, the p± spectra become harder 
and additional p±_ smearing from initial state effects becomes relatively less 
important. It is the extreme steepness of the cross sections at SPS energies that 
amplifies so greatly the sensitivity of the moderate p± yields to this aspect of 
multiple collision dynamics. 

In Fig. lb, we consider next the sensitivity of the pion yields to the sought 
after parton energy loss dynamics. The striking difference between Figs, la 
and lb is that in lb the SPS yield is not sensitive to the energy loss model in 
HIJING, while at RHIC energies the suppression of semi-hard hadrons is seen 
to be sensitive to jet quenching as predicted in Ref. ||. This seems counter 
intuitive at first because the increasing steepness with decreasing beam energy 
is naively expected to result in greater quenching for a fixed jet energy loss. 
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However, in this model the observed p± range is dominated by multiple soft 
collisions. The moderate p± < 5 GeV quarks, which fragment into the observed 
pions, are not produced in rare pQCD semi-hard interactions but are gently 
nudged several times into that p± range. Since in HIJING jet-quenching is 
restricted to only those partons that suffer a semi-hard pQCD interaction with 
p± at least pq — 2 GeV/c (the mini-jet scale), no quenching arises at this energy. 

The above conclusions are further clarified in Figs. 2 and 3, where the effects 
of p± kicks and jet quenching are shown for valence quarks and diquarks and 
gluons prior to hadronization. In Figure 2a, the valence quark plus diquark 
distribution, dN va i / dp± , is seen to be enhanced by two orders of magnitude due 
to the assumed Cronin mechanism at SPS energies whereas they are enhanced 
only by ~ 50% at RHIC energies. In Fig. 2b, on the other hand, quark jets are 
seen to be suppressed by an order of magnitude at RHIC energies, while at SPS 
energies no quenching arises because only about < 1% of the 5 GeV/c quarks 
are produced directly by hard pQCD processes. 

Figure 3a shows that the (negligible) gluon component of partons undergoing 
hard processes at SPS energies are in fact quenched by an order of magnitude 
just as at RHIC energies, but because they represent such a small fraction of the 
produced partons, their effect on the final hadron observables is negligible. Note 
the three order of magnitude rise of the mini-jet gluon component going from 
SPS to RHIC energies in Fig. 3b. In that Figure, the minimum momentum 
Pq = 2 GeV/c assumed for jet quenching is seen by the rapid drop of the gluon 
spectra beyond po- We see that the ratio of valence quarks to gluons after 
quenching at RHIC energies remains ~ 2/1 for mini-jets up to p± ~ 10 GeV/c. 

In conclusion, we found that HIJING1.35 fits the WA98 data with or without 
jet quenching. However, this model fails to account for the weak A dependence 
of the data that scale approximately with the Glauber Taa factor, and therefore 
the fit must be viewed as accidental. We note that the current WA98 data can 
also be reproduced by entirely soft physics models utilizing hydrodynamic equa- 
tions [0, , which unfortunately make no prediction for the A dependence. On 
the other hand, another Monte Carlo parton cascade model [H overpredicted 
the 7T° cross section at 4 GeV/c by a factor of 100! Those results do not prove 
the validity of hydrodynamics nor the absence of parton cascading, but empha- 
size the strong model dependence of the SPS spectra due to non-perturbative 
aspects of the problem. Those aspects, while phenomcnologically interesting, 
make it difficult to identify perturbative QCD phenomena and search for jet 
quenching. 

The main point of this work, is not to improve the HIJING soft p± phe- 
nomenology, but to emphasize the good news that at collider energies there is 
much less sensitivity to the above uncertain element of the reaction mechanism. 
At y/s > 100 AGeV the expected [[| jet quenching signature should be readily 
observable in the p± ~ 10 GeV range. Such experiments will commence in 
1999 at RHIC and should provide important tests of the theory of non-Abelian 
multiple collisions and energy loss ^ [l(| . 
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Figure 1: Invariant A + A — ► ir° cross section for central collision at SPS and 
RHIC energies are compared, a) The WA80 S + S data Q] (triangles) and the 
preliminary WA98 Pb + Pb data § (dots) are compared to HIJING1.35 
with soft p± kicks (full lines) and without px kicks (dot-dashed curves). The 
later scale with the wounded projectile number times a aa times the invariant 
distribution calculated for pp. The parton model curve from Ref. Q] is labeled 
by 'Wang'. The filled squares show pp — > tt + data scaled by the (Glauber) 
number of binary collisions times a aa for both SS and PbPb. b) Jet quenching, 
predicted at RHIC energies || , is not significant at SPS energies in the HIJING 
model. 
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Figure 2: a) Effect of soft p± kicks (g|) in HfJING on the valence quark and 
diquark dN/dp± is shown for central Pb+Pb collision at SPS versus Au+Au at 
RHIC. The dashed lines correspond to neglecting the pr kicks, b) The effect 
of jet quenching at the valence quark-diquark level with (dashed) dE/dx = 1 
GeV/fm is compared for the same systems. 
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Figure 3: a) Effect of jet quenching on dN g /dp± of mini-jet gluons (dot- 
dashed) is shown for central Pb+Pb at y/s — 17 AGeV. The solid and dotted 
curves show the dominant valence quark and diquark distribution from Fig. 2b. 
b) The effect of jet quenching in Au + Au at y/s = 200 AGeV is shown. 



10 



